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Abstract 

The empirical pairing gaps derived from four different odd-even mass staggering formulas are compared. By performing 
single-j shell and multi-shell seniority model calculations as well as by using the standard HFB approach with Skyrme 
force we show that the simplest three-point formula A^^\n) = i [B{N, Z) + B{N — 2 ,Z) — 2 B{N — 1, Z)] can provide 
a good measure of the neutron pairing gap in even-TV nuclei. It removes to a large extent the contribution from the 
nuclear mean field as well as contributions from shell structure details. It is also less contaminated by the Wigner effect 
for nuclei around N = Z. We also show that the strength of Aq (N) can serve as a good indication of the two-particle 

/o\ 

spatial correlation in the nucleus of concern and that the weakening of Aq {N) in some neutron-rich nuclei indicates 
that the di-neutron correlation itself is weak in these nuclei. 


The occurrence of a systematic odd-even staggering 
(OES) of the nuclear binding energy has long been identi¬ 
fied in nuclear physics, which is associated with the pair¬ 
ing correlation [H, 0 . It plays an important role in many 
nuclear phenomena and is the dominant many-body cor¬ 
relation beyond the nuclear mean field. Yet, in spite of 
the many efforts performed in the study of pairing corre¬ 
lations, there are still features which may be induced W 
the pairing interaction that are not well understood [3|, ^ . 
In particular, this is the case in neutron-rich nuclei, where 
the study of effects induced by pairing may shed light on 
the understanding of various exotic phenomena (see, e.g., 
Refs. dliSi). 

The simplest expression one can use to extract the em¬ 
pirical pairing gap from the OES of the binding energy is 
the three-point formula [2, |9|, which for systems with even 
neutrons acquires the form ^ 

A(3 )(]v) = -i [B{N - 1 ,Z) + B{N + 1 ,Z)- 2 B{N, Z)] 

= -^[Sn{N + l,Z)-Sr^{N,Z)] (1) 

where B is the (positive) binding energy and Sn is the 
one-neutron separation energy. The proton pairing gap 
can be defined in a similar way. The above formula indi¬ 
cates that A^^'>{N) measures the additional binding gain 
by the last neutron in the even-A system relative to the 
odd system with one more neutron. However, besides pair¬ 
ing, a number of other mechanisms may contribute to the 
OES [^, [13, 11 , I3- This includes effects induced by the 
mean field in deformed nuclei (or the Kramers degeneracy) 


and the contribution from the diagonal interaction matrix 
elements of the two-body force. 

As discussed in detail in Refs. BUI by Satula and co¬ 
workers, the contribution from the quickly varying single¬ 
particle structure of the mean field to the empirical pairing 
gap is minimized in odd-mass systems. In even systems 
where the last neutrons occupy different orbitals the single- 
particle energy contributes substantially to A*^^)(A)[3, 13|. 

Alternatively, there is another version of the three- 
point formula written as 


= ^[Sn{N,Z)-SniN-l, Z)] 

= i [B(N, Z) + B(N -2,Z)- 2B{N - 1, Z)] 

= \[S2n{N,Z)-2Sn{N-l,Z)], (2) 

which actually corresponds to A*^^) for the case of odd 
nuclei a. A(3)(jv — I) is smaller than A^^^N) in most 
cases. It is often stated that A^^'>{N — I) measures the 
pairing effect in the odd nuclei, as illustrated in Ref. Q, 
whereas A^^^N) is impacted by single particle states (see 


I, 


e.g., recent discussion in Ref. [1^ and references therein). 
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The challenging problems thus arise include: One lacks a 
measure of the pairing gap in the even system, which is 
mixed with the deformation effect, and the empirical pair¬ 
ing gap for the odd system is related to pairing only in an 
indirect manner (i.e., the energy loss due to the absence 
of pairing). The physics becomes even more obscure when 
abrupt changes occur, e.g., around shell closures. In prac¬ 
tice, as a compromise, the value of A^^\N — 1) has often 
been compared to the theoretical pairing gap calculated 
for the even systems and to the OES derived 
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from theoretical binding energies 17, [US- The direct 
comparison between the theoretical pairing gap and em¬ 
pirical OES is convenient from a computational point of 
view since only one single calculation is required, which 
avoids the complicated handling of the blocking effect in 
the odd nuclei. 

It is known that the theoretical pairing gap, e.g., that 
from the BCS theory, is not an observable and can not 
be compared with the empirical OES in a straightforward 
way. However, they can be quantitatively quite close to 
each other in most cases and both of them are still impor¬ 
tant quantities that deserves further attention. In partic¬ 
ular, within the BCS approach, the corresponding pairing 
gap is given by 

A = G'^UiVi, (3) 


where G is the pairing strength, and Ui, Vi are the stan¬ 
dard occupation numbers. This implies that the pairing 
gaps can serve as a signature of the change in two-particle 
spatial correlation/clusterization, since they are also pro¬ 
portional to '^f.UkVk- This feature is also responsible for 
the clustering of the four nucleons that eventually consti¬ 
tute the a-particle at the nuclear surface of heavy nuclei 


In this paper we would like to argue that, among the 
expressions for the OES studied here, the simple three- 
point formula A^'^(7V) removes to a large extent the con¬ 
tribution from the varying part of the nuclear mean field 
as well as contributions from other shell structure details 
and can serve as a reliable indication for the pairing effect 
in even-A systems. In other words, Aq {N) can be com¬ 
pared to the theoretical pairing gap in a semi-quantitive 
manner and thus contains fruitful information on the pair¬ 
ing effects. In particular, the abrupt changes in A^^\n) 

( 3 ) 

do have physical meaning. Moreover, by using (A) 
one can make it more convenient to extract the neutron- 
proton interaction from binding energy differences [^ . 
We will also show that A^^\n) is free from the Wigner 
effect for nuclei around N = Z and that its weakening in 
some neutron-rich nuclei may indicate that the di-neutron 
spatial correlation itself is weak in those nuclei. 

There are other formulas available for the pairing gap 
including the so-called four-point and the five-point for¬ 
mulas. The four-point formula is defined as 0 


A^^\n) = ^[-B{N + 1,Z) + 3B{N, Z) 

- 3B{N -1,Z) + B{N - 2, Z)] 

= i[A(3)(iV) + A[?)(iV)]. (4) 

That is, it measures the average value of A*^^^ in adjacent 
even and odd systems. The five-point formula is given 


by 


A(5)(A) = -[B{N + 2,Z)- AB{N + 1, Z) 

8 

+ 6B{N, Z) - AB{N - 1, Z) + B{N - 2, Z)] 

= ^[A''^\n + 2) + 2A^^\N) + aP (N)]. (5) 

The five-point formula is also used in Refs. [ii,[2i,|23. In 
Refs. @,[ 2 ^, the experimental pairing gap is taken as the 
average of adjacent ones deduced through the three-point 
formula as 


= - 


aP {N) + aP {N + 2) 


( 6 ) 


which is actually also a five-point formula involving the 
same group of nuclei as A^^) (TV) but with different weights 
for each nucleus. Our calculations show that there is no 
significant difference between the results derived from A^^e (N) 
and Ap (iV) for open-shell nuclei where the pairing gap is 
a smooth function of N. For the same reason, A'^'^\N) 
and A*^^) (A) show quite a similar behavior for most nu¬ 
clei. Noticeable differences between APe{N) and aP(N) 
may be seen where abrupt changes in pairing correlations 
are expected to happen, e.g., around shell closures, which 
is smoothed out in the former case. 
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Figure 1: Empirical neutron pairing gaps (in MeV) from different 
OES formulas for all even-even nuclei as a function of neutron num¬ 
ber N. Results which errors larger than 100 keV are marked in red 
color. The dashed curves are determined by fitting to the data. The 
solid horizontal lines show the average value of (1.08ih0.40 

MeV) and the corresponding Icr uncertainty). 


The empirical pairing gaps obtained by various OES 
formulas are shown in Fig. [1] for which the nuclear bind- 
ing energies are extracted from Refs. [2^|3^. The results 
are fitted by the expression a A“ where a and a are param¬ 
eters to be determined. Our calculations show that this ex¬ 
pression performs equally well as the usual H-dependence 
fit, as can be seen in Table [T] where the results of fitting a 
power function to the different mass formulas as a function 
of neutron number, A, and mass numbers, are given. 
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Table 1: The fittings of empirical neutron pairing gaps in even-even and even-TV odd-Z nuclei from different OES formulas as functions of A 
and TV and the corresponding root-mean-square deviation, coefficients (with 95% confidence bounds) 


Formula 

even-even 

even-even 

even-odd 

even-odd 

A(3) 

(13.43 ± 1.38)A-“'4«±u ua 

cr = 0.52 

(10.62 ± 0.8)V-‘^''^^±^ ‘^4 

cr = 0.47 

(7.58 ± 0.89)A”‘J-4i±u-‘Ja 

cr = 0.46 

(6.73 ± 0.57) 

cr = 0.42 

So 

<1 

(4.03 ± 0.46)A-'J'^’^±‘J 'J^ 

cr = 0.31 

(3.66 ± 0.33) 

cr = 0.29 

(1.40±0.27)A-'Ji^±‘J'J4 

cr = 0.29 

(1.37 ± 0.23) 

cr = 0.28 

A(4) 

(7.37 ± 0.75)A-'J'^*^±‘^'U'' 
cr = 0.35 

(6.25 ± 0.49)V-'J'^y±‘^'U'' 
cr = 0.32 

(3.46 ± 0.47)A-'J ''*^±‘^''J^ 
cr = 0.29 

(3.23 ± 0.35) 

cr = 0.28 

A(5) 

(6.93±0.70)A-u-^7±u.u2 

cr = 0.34 

(5.90 ± 0.45)V-'J'^«±‘^''J4 

cr = 0.31 

(3.94 ± 0.50)A-'J-^i±‘^''J^ 

cr = 0.26 

(3.46 ± 0.35) 

cr = 0.25 


It is also seen that the dependence of the gap npon the 
neutron or mass number is weakest for A^, (N). This is 
consistent with the conclusions of Refs, ul 3 


pairing gap should not show significant mass dependence. 


Only in nine cases one has larger than 2 MeV. 

They correspond to ®Be (4.11), ^°Be (2.57), (3.53), 

(2.8), (3.15), 20Ne (2.6), 22Mg (2.3), ^esi (2.02) 

and '^‘^Ti ( 2 . 0 ) (within parenthesis are the corresponding 
pairing gaps in MeV). 

The mean values of the pairing gaps corresponding to 
the different expressions given above are (in MeV) (N) 
= 1.46, a[?^(V) = 1.08, A(4)(Ar) = i.26 and = 


1.26. We also evaluated the empirical pairing gaps for the 
proton for which one obtains a[?^(Z) = 1.01 ± 0.36 MeV, 
which contains a smooth contribution from the Coulomb 
field. 

In Table [T] we also included calculations for the even-V 
odd-Z nuclei. However, it should be mentioned that the 


empirical OES Aq'{N) thus extracted contains a sizable 
negative contribution from the residual neutron-proton in¬ 
teraction between the odd particles in the intermediate 
odd-odd nuclei. 

As indicated in Eqs. (1), (3) and (4), all these OES 
formulas contain constantly a contribution from the mean 
field which Maks at the shell closure and persists in open- 
shell nuclei [9j . For the same reason, the OES formula Eq. 
( 2 ) may not be applicable to evaluate the pairing gaps in 
odd-V nuclei. 

The differences between the various gap formulas and 
Ap^ are plotted in Figs. [2]and[3] The dispersal of the data 
below V ^ 30 (A ^ 50) as well as around shell closures 
can be an indication of the significant mean-field contri¬ 
bution to the gaps in above regions, which is expected to 
show a A~^ dependence ll|. In this context it is worth¬ 


while to point out that the differences between A^^'>{N) 


and A[?^(V) was found to be a consequence of the gap 
between the single-particle energies of the corresponding 
neighboring orbitals [^. As can be seen in Fig. |31 the 


differences between a[?^ and A*^^) show quite small val¬ 
ues for nuclei that can be reasonably described within a 
single-j shell or nearly-degenerate systems, e.g., for nuclei 
with N between 20 and 28 {f'r/ 2 ) as well as N below ( 39 / 2 ) 



Figure 2: Differences between the neutron gaps derived from different 
gap formulas in Fig. 1 with respect to A^^(TV). The right bottom 
figure shows the neutron gaps for Z = N nuclei in comparison with 
their corresponding fitted curves from Fig. 1. 



Figure 3: Same as Fig. [2]but for the difference between A^^)(TV) and 
A^^(TV) in selected isotopic chains. 
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and above (^ 5 / 2 , 57 / 2 ) N = 50. For these systems, as we 
will show below, the difference between and is 
mainly induced by the Pauli and particle blocking effects. 

(iV) show quite large values around N = Z nuclei 
with Z = 4 — 14 and 22. One may suspect that A*^^^(iV) 
is still contaminated by the Wigner effect, which refers 
to the additional binding gained in iV = Z nuclei. This is 
indeed the case for all the other three OES formulas where, 
as seen in the lower right panel of Fig. 2, the calculated 
gap values for N = Z nuclei are systematically much larger 
than the average values. On the other hand, the calculated 

/oX 

Aq {N) values for N = Z nuclei follow nicely the average 
behavior, as it was also shown in Ref. [ 2 ^, even though 
relatively large fluctuations are present. These fluctuation 
may be due to dramatic changes in the mean field when 
going from the daughter {N — 2) system to the mother 
nucleus. For instance, they may be due to quite different 
deformation properties. An illustration of this is provided 
by the nucleus ®Be, which shows a deformed two-a-particle 
cluster structure while ®Be (as well as the mirror nucleus 
®He) is spherical. This is also the case for the nuclei ^°Ne 
and 4'^Ti which are expected to gain additional binding 
due to the strong neutron-proton quadrupole correlation. 

To illustrate the origin of the OES and the difference 
between different formulas, we start with the simple se¬ 
niority model. For a system with n identical particles in 
a single j-shell, the binding energy can be solved analyti¬ 
cally in the seniority scheme. The Hamiltonian for such a 
model can be written as [H, 


The unpaired particle blocks the scattering of other pairs 
to its own level. In practice, the odd system corresponds 
to a even system with one less particle and with a lower 
degeneracy. 

uH 

5 - 

3- 

1 - ,....,*-.-.-.'.'.3 |n/2 

_4_ 4---"'!' _I_I_^_I_I_I_I_^^ 

0 I 2 3 4 5 

Figure 4: A schematic plot for the competition between the Pauli 
blocking effect and pairing blocking effect of the unpaired odd par¬ 
ticle. 

E 2 is the energy of one pair. This term is the one that 
contributes to the theoretical OES and for systems with 
even n we have 

A[?^(n) = -iF2 = i(j + i)G. (10) 

It should be emphasized that a striking feature one thus 
finds is that the two blocking terms cancel each other in 

/oN 

Aq {n) and they do not contribute to the pairing gap. A 
schematic picture is given in Fig. |T]to illustrate this point. 
However, as implied from the figure, this is not the case 
for A^^) (n) for which one has, 

A(3)(„) = -i£;2 + iG, (11) 


H — G ^ ^m' — GS-\.S— 

m,m' >0 


( 7 ) 


where S represents the quasi-spin operator, and and a„ 
are creation and annihilation operators. G is the strength 
of the pairing interaction. The energy of the state with 
seniority v can be written as 

E{n) = -G-^{2j +3-n-v) (8) 


4 

n(n — 1) 


G-^^^^G-i(n-z;)(j + l)G 


If one assumes u = 0 for the ground state of even-even 
system and v = 1 for that of the odd system, the expression 
above can be simplified as 


£(„) = 


L2J 


(j + 1)^1 


'n' 

.2. 

( 

n' 

.2. 

- 1 

^ G -l- 

n' 

.2. 

G + 

n' 

.2. 


(9) 

E2 


where the term G/2 is due to the superfluous contribution 
from the blocking effects. This unwanted term is the main 
origin of the differences between A^^ and A(^\ as shown 
in Fig. [21 in nuclei that can be well described within a 
single j shell. 

For non-degenerate systems the pairing collectivity man¬ 
ifests itself through the correlated contribution from many 
configurations, which is induced by the non-diagonal ma¬ 
trix elements of the pairing interaction in a shell-model 
context. For two particles in a non-degenerate system with 
a constant pairing, the energy can be evaluated through 
the well known relation, 




2ji -|- 1 


= 2 . 


( 12 ) 


2ei — E 2 

L 

The corresponding wave function amplitudes are given by 

2j + l 


V = 


2ei-E2 


(13) 


where [n/2] denotes the largest integer not exceeding n/2 
and corresponds to the total number of u = 0 pairs. The 
nonlinear term in the equation above, which is propor¬ 
tional to n(n — 1) or more exactly [^] ([^] — l), is related 
to the energy loss due to the Pauli effect. The term 
in above equation indicates the energy loss in the v = \ 
odd system due to the the particle (Pauli) blocking effect: 


where is the normalization constant. All amplitudes 
Xi contribute to the two-particle clustering with the same 
phase due to the strongly attractive nature of the pair¬ 
ing interaction. The correlation energy induced by the 
monopole pairing corresponds to the difference 

Ag) = e, - \ e 2 , (14) 
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where S denotes the lowest orbital. As the gap A increases 
the amplitude Xi becomes more dispersed, resulting in 
stronger two-particle spatial correlation. This difference, 
or more exactly A — 1/2G with the self energy removed, is 
an important measure of the two-particle spatial correla¬ 
tion at the surface, reflected in a corresponding clustering 
of the two nucleons forming the pair (see, e.g.. Ref. [HI). 
This clustering induces an increase in the strength of the 
corresponding pair-transfer reaction. 

As an example, in Fig. Owe consider a simple model 
with a set of equally spaced levels with double degener¬ 
acy. We consider 16 levels and the single-particle energy 
is taken as Si = i. The total energy is obtained by diago¬ 
nalizing the corresponding Hamiltonian matrix. It is found 
that the total energy for such a system follows closely a 
relation similar to Eq. m 


E(n) 




G + 5y^l{Sb + S) 


(15) 


where Sb is the single-particle energy of the unpaired or¬ 
bital for a odd system, 6 is the energy loss due to the 
related blocking effect. It contributes significantly to the 
OES and determines its evolutions as a function of pair 
number. C/ is a parameter that is related to the pairing 
strength and level density. One has G — 1.01 and 1.06 for 
G = 0.25 and 0.5, respectively. 

Within the simple BCS context, the separation energy 
is approximately given by S 2 n{X) ss —2A and Sn{N—l) m 
—X — A, where A is the Fermi energy. One thus gets 


aP{N) PS A. 


(16) 


But one should bear in mind, as the exact solution of the 
pairing Hamiltonian in Fig. [^indicates, that there is still a 
sizable difference between the blocking effect and the OES 
especially for small systems with only few particles. 



Figure 5: Pairing gaps for a equally spaced douby degenerate 

system with 16 levels and G = 0.25 (solid circle) and 0.5 (open cir¬ 
cle). The triangles denote the contribution from the particle blocking 
effect, 5. 

From a macroscopic point of view, there may be a 
residual contribution to A^"^ from the symmetry energy 


(expected to be ne gati ve) and other non-linear terms of 
the binding energy [1^. Our calculations with the liquid 
drop model show that the average residual contributions 
are around -60 keV and -200 keV for neutron and proton 
pairing gaps, respectively. The later case is dominated by 
a smooth contribution from the Coulomb field. The main 
argument for the usage of A^®^ instead is that the smooth 
non-linear terms in the bind ing energy can be canceled up 
to the fourth order |2^, 271. That is, it removes the 
liquid-drop contribution to the OES to the largest extent. 
However, as mentioned in Ref. [l^, the disadvantage is 
that the odd-even effects may be diminished as a result of 
the averaging over nuclei further apart. 

It should also be mentioned that there is an addi¬ 
tional contribution to the binding energy from the diago¬ 
nal matrix element with m = m' (the so-called self-energy) 
which is equal to — [n/2] G. In the ideal case, one should 
have this contribution removed and the pairing gap for 
a single-j system will be of the form A = —E 2 — G12 = 
\/2{j — \/2)G. For heavy nuclei the BCS coupling constant 
roughly takes 27/A which gives a contribution to A/,' that 
is comparable but opposite to that of the smooth contri¬ 
bution from the nonlinear terms mentioned above. In the 
other words, one may expect that these two unwanted con¬ 
tributions largely cancel each other. As a result, if one as¬ 
sume that both contributions to the binding energy have 
been taken into account by the HF configuration, for the 
neutron pairing gap one is expected to have Q 


A®(iJF) =0. 


(17) 


All the variety of OES formulas studied above have 
their advantages and disadvantages (see, earlier discus¬ 
sions in Refs. Bui, 113). We will focus on the simple 
A/,' particular for two reasons: It is the one that may 
be least contaminated by the nuclear shell effect and it in¬ 
volves only three nuclei. These are important for our study 
of unstable nuclei where the experimental data are scarce 
wheras the unknown shell evolution can play a decisive 
role. For open shell nuclei, the uncertainty as related to 
the different definiation of the OES is much smaller as com¬ 
pared to the uncertainty induced by our limitted under¬ 
standing of the density functional and the shell structure 
of dripline nuclei. Our aims are to understand the pairing 
properties of neutron dripline nuclei and to extract reli¬ 
able information from binding energies based on the gap 
Aq\n). We hope that this may shed light on our un¬ 
derstanding of the stability mechanisms as well as probing 
the importance of di-neutron spatial correlations in these 
nuclei. It is theoretically suggested that di-neutron cor¬ 
relations may be enhanced in some situations such as in 
a low density region of nuclear matter and in the surface 
of finite neutron-rich nuclei [HI [33 ■ To explore this 

point further we have done systematic calculations in semi¬ 
magic neutron-rich (dripline) nuclei by using the Hartree- 
Fock-Bogoliubov (HFB) formalism, which is expected to 
be more reliable than the simple BCS approach. In the 
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standard HFB formalism, the Hamiltonian is reduced to 
the mean field in the particle-hole channel and the pairing 
field in the particle-particle channel. The HFB equations 
are, 


(H-X) A 
-A* -{H-xy 


= Ek 


(18) 


where Uk and 14 are the two components of the single 
quasi-particle wave functions. We evaluated the coordinate- 
space solution s by using the HFB solver HFBRAD in a 
spherical box 33|. In the particle-hole channel we used 


the Skyrme functional with the SLy4 parameter set 341. 
In the particle-particle channel we used the zero-range 5 
interaction given by 


Vpair{r,r') = Vo{l-r] 


P(i~) 

Po 


5(r — r') 


(19) 


Here Vq is the pairing strength, p{r) is the isoscalar nucle¬ 
onic density and po = 0.16/m“^. ij = 0 and 1 correspond 
to the volume and surface pairings, respectively. We have 
done calculations by using different pairing interactions 
but for simplicity only results corresponding to the mixed 
pairing with 77 = 0.5 are shown. The pairing strength is 
fitted to give a mean neutron gap of l.SlMeV in ^^°Sn and 
no A or isospin dependence is considered. Calculations are 
done in a box with the size r = 30 fm. Only quasiparti¬ 
cle states with energy lower than 60 MeV are taken into 
account. 

Besides OES from calculated binding energies, two dif¬ 
ferent theoretical gaps will be compared with the expert- 
mental pairing gap A^k The canonical gap A^cs, which 
consists of the diagonal elements of the pairing-field ma¬ 
trix for the lowest canonical state (LCS), and the average 
gap Amean, that is the average value of the pairing fields 
331 . 

As typical examples, in Fig. [5] we plotted Alcs and 
Amean Calculated for calcium and tin isotopes, which have 
been intensively studied quite recently from different per¬ 
spectives. The results are compared with experimental 
Aq'^ and calculated A^^ which is the pairing gap calcu¬ 
lated from the HFB binding energies by using Eq. (2). 
In the latter case, the binding energies are calculated us¬ 
ing the code HFBRAD (in the coordinate representation 
within the spherical symmetry) and compared with 
those from the code HFBTHO (in the axial deformed har¬ 
monic oscillator basis) as a cross-check for our calcula¬ 
tions. The two calculations should give very similar results 
if the influence of both the continuum and deformation is 
negligible. Indeed, as can be seen from the figure, the 
results from the two codes are very close to each other 
in most cases and come even closer if spherical symme¬ 
try constraint is applied in the latter case. In both codes 
calculations for odd-A isotopes are done within the block- 
ing app roximation using the equal filling approximation 


36l . l37l | . In practice, we have done calculations by blocking 


all possible quasiparticle orbitals around the Fermi surface 



N 


Figure 6: Calculated A^cs Amean for Ca (upper) and Sn 
(lower) together with the experimental and calculated by 

using the HFRABD and HFBTHO codes, (explained in the text). 
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and the one that gives the highest binding energy is cho¬ 
sen. A self-consistent treatment of the blocking for the 
one-quasiparticle HFB state is done in Ref. [3^ by taking 


into account beyond mean field effects using the Genera¬ 
tor Coordinate Method. In general the values of calculated 
Ap' are closer to the experimental than the two the¬ 
oretical HFB pairing gaps, which vanish at closed shell. 

The results for A^cs and Amean are practically the 
same in most cases in Fig. [B1 But it should be mentioned 
that both definitions of the theoretical pairing gap have 
their advantages and disadvantages. The use of Arcs as a 
measure of the pairing gap is questioned in Ref. [2^ and in 
Refs. [1,[I3| for dripline nuclei. In the relativistic Hartree- 
Bogoliubov calculations for selected isotopic and isotonic 
chains with a finite-range pairing force as presented in Ref. 
[T^ , it is seen that Amean are systematically smaller than 
Alcs and come closer to the A*^^) and A^^'> indicators. In 
Ref. it is also shown that the neutron A^cs vanishes 
at the drip line for nearly all semi-magic isotopic chains 
studied whereas the gap Amean can persist in some cases. 
On the other hand, Amean and A^cs can give quite similar 
prediction for bound nuclei. Calculations in Ref. Q show 
that Alcs can be closer to the experimental average gap 
Eq. (1^ t han Amean ■ In the relativistic mean field calcula¬ 
tions [17| , the calculated mean gaps and OES agree pretty 
well with each other. In Ref. [39|, the mean gap and OES 
from calculated binding energies are compared with the 
empirical OES for even-even hafnium isotopes, where it is 
found that quantatively the mean gap reproduces better 
the trend in experimental data. 

The level spacing is crucial for pairin g co rrelation which 
strongly depends on the shell structure and on the de¬ 
formation [9|. A general calculations including the defor¬ 
mation effects are done recently in Ref. with a sepa¬ 
rable pairing interaction of finite range. It is shown that 
the differences in the underlying single-particle structure 
represent the major source of uncertainty in the prediction 
of drip line and the pairing properties in neutron rich nu¬ 
clei depend substantially on the different covariant energy 
density functionals. In particular, the emergence of de¬ 
formation driving intruder orbitals can lower the chemical 
potential and make the system bound whereas extruder or¬ 
bitals with high B values can make the system unbound. A 
systematic calculation on the possible effect of the density 
dependence of the zero-range pairing force on the pairing 
correlation of dripline nuclei is also shown in Ref. . 

To explore the dependence of our calculations on the 
functionals. In Eig. [7] we extended our calculations for Ca 
and Sn isotopes shown in Fig. [5] to the neutron dripline. 
Calculations are done with the SIB Idpl . SLy4 as well 
as the recently proposed SV-bas [43| Skyrme function¬ 
als. The two-neutron dripline occurs around N = 48 and 
N = 126 for those two isotopic chains. It is thus found 
that the calculated OES agrees well with the two pairing 
gaps for known nuclei but noticable differences are seen 
as the neutorn number increases. In particular, big differ- 
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Figure 7: Calculated Alcs a-nd Amean for Ca (upper) and Sn 
(lower) isotopes upto the neutron dripline with three different 
Skyrme parameter sets in comparison with the corresponding cal¬ 
culated OES A^iN). 
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ences among different calculations are seen for Ca isotopes 
around iV = 32 in relation to the different predictions of 
the emergence of the N = 32, 34 and 40 subshells. Sig¬ 
nificant deviations between Amean and A^cs are seen in 
SV-bas calculations for nuclei beyond the dripline. In fact, 
Ref. @ has shown that there is a distinct difference be¬ 
tween the mean gap and the lowest canonical gap beyond 
drip-line for calculations with surface-peaked pairing in¬ 
teraction. A similar deviation may happen in calculations 
with the mixed pairing interaction, as in the case shown 
in Fig. [T] However, this is usually not the case for volume 
interactions. 

It should be mentioned that the time-odd field is not 
considered in the HFBTHO and HFBRAD calculations 
shown in the figure, which may affect the binding energy of 
the odd-A nucleus and the corresponding OES. The time- 
odd field can contribute to the ground state mean fields 
of odd-A and odd-odd nuclei (breaking the Kramer’s de¬ 
generacy) as well as in nuclear rotation and other dynamic 
processes [i^. The time-odd field, which is still relatively 
poorly understood and functional dependent, is neglected 
in most relativistic and non-relativistic mean field calcu¬ 
lations. The possible influence of the time-odd field in 
light N = Z nuclei is discussed in Ref. [i^. System¬ 
atic Skyrme Hartree-Fock plus BCS calculations with the 
time-odd field for odd-A nuclei betweenlG < Z < 92 are 
presented recently in Ref. [d^. It was shown that the 
influence of the time-odd field is generally small and de¬ 
crease rapidly with increasing mass number. In particular, 
there is no indication that the deviation between experi¬ 
mental and theoretical OES can be improved by the inclu¬ 
sion of the time-odd field [d^, from which the fluctuation 
induced is much smaller than that from the different func¬ 
tional. Skyrme HFB calculations with the time-odd field 
are also done for the whole mass table in Ref. and for 
rare earth nuclei with 63 < Z < 75 and 78 < N < 104 
in Ref. [s^. It is found that the average energy shift as 
induced by the time-odd field is only 50 keV but with a 
large deviation of 42 keV and hence the equal filling ap¬ 
proximation is precise for most practical calculations. A 
similar conclusion may be draw from Gogny HFB calcula¬ 
tions where full calculations with the time-odd field give 
very similar results to those from the filling approximation 
and the inclusion of the time-odd field hardly improves the 
gap description [i^. Recent systemtic calculations on the 
effect of the time-odd field in the relativistic mean field 
approach can be found in Refs. 4^, 5^, 51| and references 
therein. Ref. [d^ found that the time-odd field in the 
relativisitic mean field approach always induces an addi¬ 
tional binding but weakly affects the relative energies of 
different quasiparticle states in medium and heavy mass 
nuclei. The influence of that attractive contribution from 
the time-odd mean field to the odd nuclei into OES can be 
of up to 10% in light nuclei and of around 5-6% in heavy 
nuclei. However, even in that case, the required modifica¬ 
tion of the strength of pairing force is modest [d^, which 
is significantly smaller than earlier expectations. In Ref. 



R [fm] R [fm] 


Figure 8: Upper: The two-neutron spatial correlation plots for ^®Ca 
(left) and ®‘^Ca (right). Lower: Same as upper but for ^^®Sn (left, 4 
holes) and ^®®Sn (right, 4 particles). Calculations are done with the 
SLy4 force. Notice that the scale is different. 


[ 2 ^ . it is suggested that the Ap^(iV) can be a measure of 
the sole pairing gap if the contribution from the time-odd 
reversal symmetry breaking field cancels out the contribu¬ 
tion from the smooth part of the total energy. 

In Ref. @ it is mentioned that for spherical nuclei 
where the drip line occurs around neutron closed shells 
the pairing gap is expected to be reduced. Our calcula¬ 
tions show that this is indeed the case for neutron-rich Ca 
and Sn isotopes. Our understanding of the shell structure 
in Ca isotopes have been significantly extended recently 
(see, e.g.. Refs, [s^, [H, [H^). N = 32 and 34 have been 
shown to be new magic numbers in Ca isotopes, a con¬ 
clusion which is supported by HFB calculations . The 
A^^\n) values for ^°’^^Ca are much smaller than those 
for nuclei in the O/ 7/2 shell. The pairing gap for ^"‘Ca is 
expected to be even smaller. In the upper panel of Fig. [51 
following the same procedure as in Ref. [2l| , we evaluated 
in the canonical basis the two-neutron spatial correlations 
in ^'^Ca and compared it with that in ^®Ca. The results 
are plotted as a function of the angle 6 between the two 
neutrons and the radius R = ri = r 2 where r is the radius 
of the single-particle wave function. For strongly corre¬ 
lated wave functions, the two-neutron spatial correlation 
is expected to peak at 0 = 0 at the surface. As can be seen 
from the figure, the peak in '^^Ca is much stronger than 


that 


56 


Ca. 


The spatial correlation of the neutron pair can also be 
inferred from the pairing/abnormal density p{r). in Fig. 
O we plotted the calculated pairing tensor which is defined 
as K,{r) = r^p{r) and is related to the pair transfer form 
factor, for two nuclei ^®Ca (solid line) and ^®°Sn. ^®Ca 
is at the neutron dripline for which the calculated pairing 
tensors and pairing gaps from difference functionals can be 
quite different. The nucleus ^^'^Sn is within the dripline. 
The pairing tensors calculated from the three functional 
are quite close to each other. But as can be seen from Fig. 
[71 there is a significant difference between the theoretial 
mean gap Amean and the OES from the theoretical bind- 
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Figure 9: The pairing tensor as a function of radius r for neutron- 
rich nuclei ^®Ca (solid line) and ^®*^Sn (dashed line) as calculated 
from different Skyrme parameterizations. 


ing energies. A similar deviation is also seen in Ref. [27 1. 


The origin for such a large deviation is not exactly known 
yet. It may be related to the fact that, for Sn isotopes in 
this region, the chemical potential A is pretty close to zero 
and the quasiparticle orbitals involved in the blocking cal¬ 
culations are all unbound. For such cases, the theoretical 
OES and gaps have to be carefully interpreted. 

There is a noticeable kink at TV = 66 for tin isotopes 
which may be related to the occupancy of the low degener¬ 
acy orbital Si/ 2 - A sudden drop occurs in the experimental 
pairing gaps in the Sn isotopes when going from TV = 82 to 
N = 84. This can be reproduced by the HFB calculations 
with the volume and mixed pairing interactions Q. That 
drop is related to the reduced level density above N=82. 
That is, there is a noticeable gap between the I/7/2 orbital 
and the higher- l ying ones, which may result in a = 90 
sub shell closure [52j| . The reduced pairing collectivity in 
these nuclei can be clearly seen in the two-neutron spatial 
correlation plots given in the lower part of Fig. |8j where we 
compared ^^®Sn, which have four holes in the core ^®^Sn, 
and ^®®Sn, with four particles above the core. As can be 
seen from the figure, the peak in ^^®Sn is around one order 
of magnitude stronger than that in ^®®Sn. 

HFB describes very well the pairing gap as well as 
the two-particle spatial correlation in open-shell nuclei. A 
known problem is that both the BCS and the HFB conden¬ 
sates collapse at closed shell, since in this case all orbitals 
are fully occupied (c.f., Eq. (13)). In fact, as the sys¬ 
tematic calculations in Ref. [I 3 shows, roughly 20-30% of 
the known nuclei contain collapsed BCS or HFB conden¬ 
sates. It is therefore fair to affirm that the failure of these 
two approaches to reproduce the OES does not necessary 
mean that there is no pairing correlation in these nuclei. 
Two typical examples are "'^®Ca and ^®^Sn shown in Fig. 
m These two nuclei show values similar to neighbor¬ 
ing nuclei below the closed shell. Indeed our calculations 
with the seniority model [s^ show that the two-neutron 


spatial correlations in these two nuclei are as strong as 
those in neighboring ones. Again, is a good indica¬ 
tion of the two-neutron correlation in these cases. Anyway, 
as mentioned above, comparisons between the theoretical 
gaps from even-even nuclei can only be compared with the 
empirical pairing gap in a semi-quantitative manner. 

In summary, in this paper we compared the pairing 
gaps derived from four different OES formulas. We showed 
that A)^' gauge very well the nuclear pairing correlation 
since it removes to a large extent the contribution from 
the single-particle structure of the nuclear mean field and 
the shell effect. It can serve as a reliable filter if one is pri¬ 
marily interested in evaluating the pairing effect in struc¬ 
ture model calculations. This is particularly interesting 
for the study of dripline nuclei since experimental data 
on those nuclei are usually scarce and show large fluctua- 
tions of shell structure. In addition, AJ^^ is expected to be 
less contaminated by the Wigner effect for nuclei around 
N = Z. We have also shown that the strength of A® can 
be a good semi-quantative indication for the two-particle 
spatial correlation. This is supported by our calculations 
with the HFB model with Skyrme force as well as with 
the multi-shell seniority model. Moreover, we found that 
the weakening of [N) in some neutron-rich nuclei in¬ 
dicates that the di-neutron correlation is weak in those 
nuclei. As examples, the pairing gaps and di-neutron spa¬ 
tial correlations in neutron-rich calcium and tin isotopes 
are evaluated in detail. Calculations for the different pair¬ 
ing gaps and OES with different functional quantatively 
agree well with each other as well as with experimental 
data. However, these results need to be carefully examined 
when one approaches the dripline, where large deviations 
between different functionals and different definations of 
the pairing gaps may be seen. 
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